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ABSTRACT 


At a Prandtl number of about 0.2, the predictions of accepted 
correlations for heat transfer in fully established tube flow differ by a 
factor more than two. Since gas turbine systems proposed for propulsion 
utilize working fluids with the Prandtl number in this range, it was 
necessary to resolve the discrepancy. By mixing helium with xenon or 
hydrogen with xenon, the range 0.16 < Pr e 0.7 can be obtained. 
Measurements with these mixtures in a vertical tube showed that the Colburn 
analogy and Dittus-Boelter substantially overpredict the Nusselt number for 
constant property conditions; best agreement was provided by relations 
suggested by Petukhov and by Kays. For moderate variation of gas 
properties (1< T /t, < 2.2) the correlation for average friction factor by 
Taylor was vavitied: the exponent on the Prandtl number in his equation for 


heat transfer was modified to 0.65 in order to accommodate these new data. 
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NOMENCLATURE 


Flow area of test section; 
Velocity of sound; 


Specific heat at constant pressure; 
Specific heat at constant volume; 


Inside diameter 
Acceleration due to gravity; 


Dimensional conversion factor; 


Average mass flux, m/A; 
Heat transfer coefficient; 
Specific enthalpy; 


Thermal conductivity; 


Mass flow rate; 

Distance from start of heating to PT #2; 
Molal mass; 

Pressure; 

Heat flux; 

Gas constant for a particular gas; 
Universal gas constant; 

Temperature ; 


Axial distance from start of heating; 
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NOMENCLATURE (continued) 


Greek Symbols 


— e/K Force constant in Lennard-Jones potential; 
; Y Ratio of specific heats, e/ey3 
bad u Absolute viscosity; 
v Kinematic viscosity; 
fe) Density; 
( 
: fof Force constant in Lennard-Jones potential; 


Nondtmenstonal Parameters 


Peet eee F 


2 
f Friction factor, g. Py avo Per /2 LG ; 
G. Grashof number based on wall heat flux, 
e g p* "1 kT), 
Nu Nusselt number, h D/k ; 
f é 2 
: P Pressure drop, QP, 8. (Py - p)/G” ; 
Pr Prandtl number, cyn/k $ 
( 
: + " ; 
q Heat flux parameter, q,,/6 Cot T,) ; 


Reynolds number, G D/u ; 


Subsertpte 


av 


DB 


fr 


NOMENCLATURE (continued) 


Lengthwise average; 
Evaluated at bulk temperature; 


Constant property condition; 


Dittus~-Boelter; 
Frictional; 


Hydrogen gas; 


Inlet; an index; 
Modified; 

Reference; 

Wall; 

Xenon gas; 

At start of heating; 


At pressure tap 2. 
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INTRODUCTION 


The use of gas mixtures has been proposed for use in closed cycle 
gas turbines with applications in power systems undersea and in space as 
well as for propulsion of ships, aircraft, buses and rail units. Helium 
and hydrogen both possess excellent heat transfer properties but require 
many stages of turbomachinery in closed cycle gas turbines. The addition 
of higher molecular weight gases to helium or hydrogen to form binary 
mixtures can reduce the number of stages in turbomachinery but increases 
heat exchanger surface area because of the lower thermal conductivity of 
the heavier gases. The Prandtl number can be as low as 0.16 for these 
mixtures. In the past, most designs utilizing these gas mixtures were made 
by using correlations of the Colburn [1933] or Dittus-Boelter [1930] type 
which had resulted from tests using gases with a Prandtl number around 0.7, 
or with liquids. 

In fact, some well-known texts do not recognize that a Prandtl 
number less than 0.6 can occur for gas mixtures. Several relations have 
been reported which differ considerably from Colburn and are suggested for 
use over a wide range of Prandtl numbers, some as low as 0.1 (Kays, 1966; 
Petukhov, 1970; Sleicher and Rouse, 1975; Gnielinski, 1975; and Churchill, 
1977]. The relations recommended by these investigators are tabulated in 
Table 1. The large variations in Nusselt number resulting from these 
relations are shown in Fig. 1 as a function of Prandtl number. At a 


Prandtl number of 0.2 there is a factor of 2.2 between the lowest and 
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Figure 1. Comparison of correlation equations recommended 
by several investigators. 


highest value recommended! The danger of extrapolating heat transfer 
correlations based on data at Pr = 0.7 is emphasized by Pierce [1981] in 
an application to the closed gas turbine cycle. 

The authors of this paper have used a number of binary mixtures of 
light and heavy gases to test the correlations shown in Fig. 1; for later 
reference, equation numbers are listed on the table. Pickett, Taylor and 
McEligot [1979] employed mixtures of helium and argon in order to obtain 
Pr = 0.42 and 0.49 while Serksnis, Taylor and McEligot [1978] conducted 
measurements at Pr =~ 1/3 by using a mixture of hydrogen and carbon 
dioxide. In the present investigation, additional data have been obtained 
in the range 0.18 < Pr < 0.3 by mixing hydrogen with xenon and helium 
with xenon. Table 2 summarizes the range of data obtained in this 


investigation. 
TRANSPORT PROPERTIES OF THE MIXTURES 


The properties needed for this study were compressibility, 
viscosity, thermal conductivity, specific heat, enthalpy, speed of sound, 
and the gas constant. The properties of air have been studied extensively; 
for the comparison data the NBS tables of Hilsenrath et al. [1955] were 
used in this investigation, The properties of helium-xenon mixtures were 
calculated theoretically. For all gases the viscosity and thermal 
conductivity were assumed to be independent of pressure. 

The helium-xenon mixtures were assumed to be ideal gases, thus 
making the compressibility equal to unity. This assumption is reasonable 


for the range of pressures (260 - 1000 kPa or 2.6 - 9.9 atm) and 
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temperatures (290- 980°K or 62 - 1310°F) used in this experiment. Since 
helium and xenon are monatomic and the temperature range in this study was 


not too large, the relation [Reynolds, 1968] 


iS c, = (3/2) R= (5/2)R/M (8) 


was used to calculate the specific heat. 

Using the ideal gas and constant specific heat assumptions, one may 
derive simple equations for the enthalpy and speed of sound [Reynolds and 
Perkins, 1968] 


izc = (T-T ) 
and p ref (9) 


c= 7yRT = V7(5/3)RT 

The Lennard-Jones (6~12) potential can be employed in the Chapman-~ 
Enskog kinetic theory to predict thermal conductivity, viscosity and 
Prandtl number of binary mixtures of inert gases [Hirschfelder, Curtiss and 
Bird, 1964]. There has been considerable experimental study of the pure 
gases but, unfortunately, few data exist on the mixtures. 

With force constants, €/k and 5, suggested by Hirschfelder, Curtiss 
and Bird [1964], the predicted viscosity for pure helium falls about eight 


percent below the data of Clarke and Smith [1968], Dawe and Smith [1970] 


and Kalelkar and Kestin [1970] at temperatures around 900°C (1650°F). 
Likewise, the predicted thermal conductivity is about nine percent lower 
than the measurements of Saxena and Saxena [1968] up to 1100°C (2011°F). : 


In contrast, using the force constants suggested by DiPippo and Kestin 


[1969] leads to essential agreement with the values recommended by the 


Thermophysical Properties Research Center [Touloukian and Ho, 1970]. 

For pure xenon, force constants from Hirschfelder, Curtiss and Bird 
[1964] predicted viscosity that falls about five percent below the data of 
Dawe and Smith [1970] and Kestin, Ro and Wakeham [1972]. The predicted 
thermal conductivity is about twelve percent lower than measurements of 
Saxena and Saxena [1969] up to 1200°C (2191°F). DiPippo and Kestin did not 
report force constants for xenon, but Kestin, Ro and Wakeham suggested 
force constants that predicted values of viscosity that are in very good 
agreement with the experimental data and predicted values of 
thermal conductivity that are about four percent lower than the reported 
measurements. 

The properties of the helium-xenon mixtures are shown in Fig. 2. 
The solid curves were calculated using the force constants recommended by 
DiPippo and Kestin [1969] for helium: o = 2.158 A and e/k = 86.2°K and 
those recommended by Kestin, Ro and Wakeham [1972] for xenon: o = 2.858 A 
and €/k = 285.2°K. 

The viscosity of the mixtures varies considerably with the 
molecular weight of the mixture. At 277°C (530°F) the maximum viscosity of 
the mixture is 41% higher than that of helium and five percent higher than 
xenon. Viscosity predicted using constants from DiPippo and Kestin [1%9] 
and Kestin, Ro and Wakeham [1972] agrees to within three percent of the 


data of Trautz and Heberling [1923] and Thornton [1960]. As with pure 


gases the viscosity increases with temperature. 
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The mixture thermal conductivity decreases by a factor of 23 as the 
molecular weight increases from pure helium to pure xenon over the range of 
mixture temperature in this investigation, and mixture conductivity also 
increases with temperature. Agreement with the data of Thornton [1960], 


Gandhi and Saxena [1968] and Mason and von Ubish [1960] is good for 


temperatures ranging from 18° to 90°C (65° to 194°F). The only data 
available at higher temperatures are those of Mason and von Ubish at 520°C 
(98°F) and these are almost 10% higher than the predictions; however, in a 
critical review Gandhi and Saxena [1968] have observed that the 
measurements of Mason and von Ubish appear to be systematically higher than 
others they reviewed. 

As a consequence of the variation of thermal conductivity and 
specific heat vs. molecular weight the Prandtl number decreases to a 
minimum of about 0.21 at M ~ 50 from about 0.667 for the pure gases. It 
is about 0.23 at the molecular weight of air. 

It is interesting to note that even though the force constants of 
Hirschfelder, Curtiss and Bird [1964] did not adequately predict thermal 
conductivity and viscosity for pure helium and xenon, their predictions for 
the mixtures did not differ greatly from the ones used in this 
investigation. 

The hydrogen-xenon mixtures were also assumed to be ideal gases, 
making the compressibility equal to unity. This assumption is also 
reasonable for the range of pressures (750 - 800°KPa or 7.4 - 7.9 atm) and 
temperatures (290 - 830°K or 62 - 1040°F) used in this experiment. The 


specific heat of this mixture was calculated by the relation 


ll 


. 
ee Ge On ae ea io 


Csatcure = (mass fraction ae + (mass fraction, ,) c 


2 p,Xe 


The enthalpy and speed of sound were calculated in the same manner as with 
the helium-xenon mixture. 

The hydrogen-xenon mixture properties are shown in Fig. 3. The 
solid lines for viscosity were calculated using the viscosities of hydrogen 
and xenon and the method recommended by Hirschfelder, Curtiss and Bird 
[1964] for calculating mixtures of monatomic gases. Hydrogen is, of 
course, polyatomic and predicting its properties is more complex than 
predicting those of monatomic gases, The solid lines for thermal 
conductivity were calculated using the method of Lindsay and Bromley 
[1965]. 

The only measurements of viscosity of hydrogen-xenon mixtures found 
in the literature were those of Trautz and Heberling [1934] which are in 
good agreement with these predictions, especially at the mole fraction of 
xenon (0.21) in the mixture used in the present investigation. The 
viscosity of the mixture increases with temperature and the increase in 
viscosity from pure hydrogen to pure xenon is threefold. 

The mixture thermal conductivity decreases by a factor of 30 as the 
molecular weight increases from pure hydrogen to pure xenon over the 
temperature range of this investigation. Only two investigations of 
thermal conductivity of hydrogen-xenon mixtures have been reported. Barua 
[1960] measured conductivities for mixtures of eight volume-percentages of 


xenon from 0 to 100% at temperatures of 30°C (85°F) and 45°C (111°F). 
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Saxena and Tondon [1971] measured conductivities for five mole fractions of 
xenon from 0 to 1.0 at 40°, 65° and 93°C (104°, 149° and 199°F). The 
measurements of the two investigations are in good agreement with each 
other, particularly near the mixture concentration (mole fraction of 
Xe = 0.21) used in this investigation, but are as much as fifteen percent 
higher than predicted values at 93°C (199°F). 

For the present concentration, the experimental values of thermal 
conductivity were extrapolated to the maximum mixture temperature used, and 
it was found that the measured conductivity was twenty percent higher than 
the predicted value, A method recommended by Kestin [personal 
communication, 1983], but not with complete confidence, is described by 
Clifford et al. [1980]. The thermal conductivities calculated by this 
method are closer to the experimental values than those predicted by the 
method of Lindsay and Bromley, but are still fifteen percent lower at 
283°C. As a consequence of this discrepancy in thermal conductivity, the 
heat transfer data of this investigation were reduced twice, once using the 
lower predicted values and once using the values obtained by extrapolating 
the experimental values. The normalized Nusselt number was found to be six 
to ten percent higher when the calculated thermal conductivity was used 
rather than the extrapolated experimental value. 

As a result of the variations of the calculated thermal 
conductivity, viscosity and specific heat versus the molecular weight, the 
calculated Prandtl number decreases to a minimum of about 0.18 at M = 47 
from 0.707 for hydrogen and 0.667 for xenon. It is about 0.20 at M= 29 


which is the molecular weight of both air and the hydrogen-xenon mixture 
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used in this investigation. When based on the experimental thermal 


conductivity a value of 0.18 is predicted for this mixture. 


Apparatus 

The experimental apparatus, arrangement and procedure were similar 
to those used by Park, Taylor and McEligot [1982]. In the present 
investigation the loop was a closed circuit, as shown in Fig. 4, due to the 
extremely high cost of the xenon gas. A single-acting "Gas Booster Pump” 
from Haskel Manufacturing Company circulated the gas mixtures through two 
pressure regulators, a plenum, a cooler which removed the heat of 
compression, a tubular flowmeter, the instrumented test section, another 
cooler to remove the energy added in the heated test section, a UGC 
densitometer, another plenum, a control valve and then back into the pump. 
The two pressure regulators and plenum were installed to remove the 
Pressure fluctuations in the flow created by the pump. Two pressure 
transducers were used to measure the pressure fluctuations; a Model SCD 147 
from Data Instruments, Inc. was located just beyond the first cooler and a 
Kulite Model XT-140-100G subminiature pressure transducer was mounted flush 
with the inside of the tube immediately beyond the elbow at the entry of 
the test section to measure the pressure fluctuations of the flow entering 
the test section. 

The vertical test section was a circular tube of Inconel 600 with 


an inside diameter of 5.87 mm (0.231 in) and a wall thickness of 0.28 mm 
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(0.011 in). It consisted of a heated section 60 diameters in length 


preceded by an unheated section of 56 diameters which insured that the flow 
approached a fully developed velocity profile prior to heating. The test 
section itself served directly as an electrical resistance heater. 
Alternating current from a line voltage stabilizer was supplied via 
variable transformers to the test section through thin stainless steel 
electrodes brazed to the tube. The current was measured using a Weston 
current transformer and a Weston Model 370 ammeter. A high impedance Fluke 
differential voltmeter measured the voltage across the test section to 
serve as a check on power measurements. Three pressure taps, with holes of 
about 0.30 mm (0.012 in.) were used. One was located 5 diameters below the 
lower electrode and the other two were 50 and 54 diameters above it, near 
the upper electrode. The static pressure at the test section inlet was 
measured with a Heise bourdon tube gauge and the pressure drop between taps 
was measured with an MKS Baratron Type 77 Pressure Meter and Head. The 
fluctuating signals from the Kulite pressure transducer were recorded on a 
Hewlett-Packard x-y recorder. 

Twenty-one premium grade Chromel-Alumel thermocouples, 0.13 om 
(0.005 in.) in diameter were spot welded to the heated section of the tube 
using the parallel junction technique suggested by Moen [1960]. Thermal 
conduction error was calculated from the heat loss calibration data and a 
relation developed by Hess [1965]; the thermocouple conductance was 
estimated from the emissivity of the bare wire and a natural convection 
correlation for small Rayleigh numbers [Kreith, 1973]. The correction was 


of the order of 1% of the difference between the tube temperature and the 
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environmental temperature; at Re~ 3 x 104 this correction was equivalent 


to 1-1/2 to 2% in the Nusselt number. 


Procedures and Preliminary Measurements 

The test section was a bare tube surrounded by a draft shield so 
the heat loss was by radiation and natural convection. The heat loss was 
calibrated as a function of axial position and temperature from 
measurements with a vacuum on the inside of the test section. the 
uncertainty of the heat loss measurements was estimated to be about 1%. 
The heat loss data were represented by cubic equations that deviated no 
more than three percent except at very low test section temperatures 121°C 
(250°F) where the deviation might be as much ds ten percent. For tests 
with gas flow the heat loss was usually less than 10% of the heat 
generated. A few runs with inlet Reynolds number less than 36000 and 
maximum wall temperatures greater than 650° (1200°F) had heat losses 
exceeding 20% of the heat generated. The worst case was a heat loss of 412% 
for helium-xenon (M = 83.8) with an inlet Reynolds number of 34000 and a 
maximum wall temperature of 680°C (1256°F). The electrical resistance of 
the test section was also measured as a function of temperature during the 
heat loss runs; the estimated uncertainty was also about 12%. 

The closed loop was pressurized with gases and gas mixtures from 
high pressure cylinders. A cylinder of helium-xenon gas mixture with a 
molecular weight of 83.8 was obtained from the NASA-Lewis Research Center. 
Helium-xenon mixtures of lesser molecular weight were obtained in the loop 


by adding helium to the existing gas mixture. On the other hand, the 
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hydrogen-xenon mixture was obtained by mixing high purity hydrogen directly 
with high purity xenon. The partial pressures of the gases in a mixture 
were used to calculate the quantity of helium needed to reduce the 
molecular weight of the helium-xenon mixture and also to determine the 
amounts of hydrogen and xenon needed to obtain the particular mixture for a 
given run. 

The UGC gas densitometer was calibrated with air, argon, helium and 
the helium-xenon gas mixture (M = 83.8) over the pressure range of this 
investigation. The densitometer was used to verify the procedures for 
mixing the gases to predetermined concentrations in situ and also for 
checking the concentrations after a series of runs. This was done by 
comparing the measured density with the density calculated using the 
measured temperature and pressure and the perfect gas law. The measured 
and calculated densities usually agreed within one percent and the mixture 
molecular weight did not change measurably during the runs. Mass flow rate 
was determined with the tubular flowmeter, which was itself calibrated over 
the range of interest using several positive displacement meters in 
parallel. The mass flow rate could be determined within an uncertainty of 
1.5% or better. 

Park, Taylor and McEligot [1982] reported the results of an 
experiment on the effects of pulsating flow on heat transfer with air with 
the same apparatus modified to an open loop configuration, so that the mass 
flow rate could be measured directly with positive displacement meters at 
the exit of the loop. The static pressure fluctuations were as large as 


could be attained with the gas booster pump in this system (9 to 35%). 
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Frequency ranged from 2.1 to 3.6 Hertz, q* from zero to 0.0034 and Reynolds 
Numbers from 18,000 to 102,000. For the range of conditions in the present 
investigation the effect of pressure fluctuations from 26 to 35% had less 
than a 2% effect on the Nusselt number [Park, Taylor and McEligot, 1982]. 
The pressure fluctuations in the gas mixtures were never more than 0.62% 
(usually less than 0.22) and had a frequency no greater .° Wn 2.0 Hertz. 
Thus, any effect of pressure fluctuations on the Nusselt number of the 
present gas mixture experiments is believed to be negligible. Results from 
runs made with air pulsating at the present levels in the closed loop 
compared closely with those from similar runs conducted with the loop open 
and no pressure fluctuations, 

The procedure for the experimental runs was to introduce the proper 
amounts of the gases into the closed loop and then circulate the mixture 
until the density of the gas reached steady state as measured by the 
densitometer. The mass flow rate was set to give the required Reynolds 
number, and then the electrical power to the test section was adjusted to 
give a series of maximum wall temperatures of approximately 120°C (250°F), 
260°C (500°F), 400°C (750°F), 540°C (1000°F) and 680°C (1250F). After 
covering the range of wall temperature the mass flow rate could be changed 
to give the next required Reynolds number. Once the range of Reynolds 
number was covered, the power was shut down and the gas flow stopped. At 
this time the density was again measured to determine whether the molecular 
weight had changed due to preferential leakage of the hydrogen or helium. 
No change in molecular weight could be detected at any time during the 


experiment. 
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During each run the wall temperatures along the heated and unheated 
test section were recorded along with the inlet static pressure and the 
pressure difference across both the tubular flow meter and the heated test 
section. The current through the test section and the voltage drop were 
recorded. The bulk temperature of the gas entering the tubular flow meter 
was measured and the temperature of the gas entering the heated test 
section was deduced from this measured temperature and the wall 
temperatures of the unheated section just downstream of the start of 
heating. 

The data were reduced to give local heat transfer and fluid flow 
parameters. As described later in the Experimental Results, the ratio of 
measured Nusselt number to the Nusselt number predicted by Dittus and 
Boelter [1930] was plotted as a function of the ratio of wall-to-bulk gas 
temperature and was then extrapolated to T,/Tp = 0 to deduce a constant 
property Nusselt number versus position for each gas mixture and Reynolds 
number. This constant property Nusselt number was compared to the 
predicted values suggested by several investigators mentioned earlier. The 
heat transfer and friction data with property variation were reduced and 


compared with existing correlations. 


Experimental Uncertainties: 

The experimental uncertainties were estimated by the method of 
Kline and McClintock [1953]. Typical uncertainties for the Nusselt number 
were about 8% at x/D ~ 1.3 decreasing to 5% at x/D > 24 for low heating 


rates, and about 1.4% at x/D ~ 1.3 increasing to 4% at x/D < 24 for the 
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higher heating rates. These estimates are in good agreement with the 
estimates made by Serksnis [1977] for Hy - CO, experiments and by Pickett 
[1976] for He-Ar experiments, both in open loop configurations. For the 
low heating rates, the dominant uncertainties are provided by the inlet 
bulk gas temperature and the wall temperatures. For higher heating rates, 
the uncertainty increases with x/D because the uncertainties in tube wall 
temperature and, therefore, temperature difference increase significantly 
with temperature level, while the contributions of uncertainties in mass 
flow rate, electrical power and inlet gas temperature remain small. 
Typical uncertainties for air, He-Xe M ~ 40) and Ho-Xe (M 229) data are 


shown in Table A2 in Appendix A. 


Reproducibility 

The reproducibility of the measurement technique was checked in two 
ways. Air data in steady flow had been obtained previously in two other 
test sections by Serksnis, Taylor and McEligot [1978] and Pickett, Taylor 
and McEligot [1979]. It was found that each had a series of runs at Re; 
near 80,000 and various heating rates, so these were compared to present 
measurements at the same conditions. For the three sets of data which 
spanned a five-year period, it was found that in the downstream region the 


0.8 pp0.4) agreed to within three 


normalized Nusselt number, Nu/(0.021Re 
percent at low heating rates (T/T, ~ 1.2) and within two percent at higher 
heating rates (1.4 < T/T, < 1.8). 


Secondly, the reproducibility of the present measurements was 


tested at the end of the experiments by duplicating one of the first runs 
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with Re; ~ 60,000, q* = 0.0014 (maximum T,,/T, ~ 1.5) and steady conditions. 
| The mass flow rate could be reproduced to better than 0.2%, the test 
section inlet pressure to within less than 0.12 and the electrical current 
within the accuracy of ammeter (~0.25%), The resulting values of 


, oe (TY - Th,in)max differed by 2.1%, leading to agreement of the fully 


developed Nusselt numbers within less than 3% again. 


{ EXPERIMEETAL RESULTS 


The ranges of the present data have been previously presented as 
_ Table 2. Tabulations of these data are provided in Appendix B. Data for 
helium~argon mixtures are tabulated in the report by Pickett [1976] and for 


hydrogen-carbon dioxide mixtures by Serksnis [1977]. 


Heat Transfer with Constant Properties 

To compare experimental data with constant property correlations 
the measurements were extrapolated to the constant property idealization by 
an approach like that of Malina and Sparrow [1964]. For this method, a 
series of experimental runs was taken with the same inlet Reynolds number 
and gas composition but with successively higher heating rates. At each 
thermocouple location the normalized value of the measured Nusselt number, 
Nu/Nup,, was plotted versus the local ratio of wall to bulk temperature, 
T,/Tp. An extrapolation to T,,/T, = l yielded the deduced Nusselt number 


for constant property conditions, Nu, 


p* 


For fully developed conditions with air as the fluid, the ratio of 
Nucp/Nupg was about 0.94 for Reynolds numbers from about 34,000 to 85,000 
without any apparent dependence on Reynolds number. The variations of 
Nucp/Nupg with T/T, for several thermocouple locations are shown in 
Figure 5 for air and three binary gas mixtures. It can readily be seen 
that the extrapolated values of Nu,p/Nupg decrease from 0.94 for a Prandtl 
number of 0.72 to 0.86 at Pr = 0.42, 0.76 at Pr = 0.21 and 0.66 at the 
lowest Prandtl number. This reduction in the normalized value corresponds 
to the trend predicted by an analysis by Pickett, Taylor and McEligot 
[1979]. 

The constant property Nusselt numbers for Prandtl numbers from 0.18 
to 0.72 for three Reynolds numbers are compared to various correlations 
from references from Dittus-Boelter [1930] through Churchill [1977] in 
Figure 6. The solid symbols denote data from the present investigation. 
The air data from Pickett [1976] and Serksnis [1977] and the present 
investigation agreed to within less than 32 and are represented by a single 
solid symbol. As the Reynolds number increases, tthe deviation between the 
predictions of various correlations also increases. At Reynolds numbers of 
34,000 and 60,000 there is little difference between the values predicted 
by Kays [1966] and Petukhov and Popov [Petukhov, 1970], and both are in 
reasonably good agreement with the experimental data for Prandtl numbers 
from about 0.2 to 0.72. However, for a Reynolds number of 84,000, the 
difference between th Petukhov and Popov and Kays’ correlations widens and 


the Petukhov relation agrees more closely with experiments at the lower 


Prandtl number. 
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Figure 6a. Comparison between measured Nusselt number and Nusselt 
number predicted by correlations proposed by other 
investigators for constant properties. 
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Figure 6b. Comparison between measured Nusselt number and Nusselt 
number predicted by correlations proposed by other 
investigators for constant properties. 
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Figure 6c. Comparison between measured Nusselt number and Nusselt 
number predicted by correlations proposed by other 
investigators for constant properties. 
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There could be concern about the possibility of thermo-diffusion at 
high temperatures with the high ratios of molecular weights of the gases 
used in the binary mixtures. However, any effect of thermo-diffusion 


should reduce to zero as the data are extrapolated to a T/T, = 1. 


Heating with Property Variation 

If power densities are to be increased, or the weight of systems 
such as the closed cycle gas turbine reduced, moderately high heating rates 
must be employed. Then, the temperature-dependence of the fluid transport 
properties will cause significant variation in properties appearing in the 
correlation equations for both heat transfer and friction coefficients. 
The results and correlations based on the constant properties idealizations 
become invalid. In this section the proposed modifications accounting for 
property variations are examined. 

Since pressure taps were installed only near the entrance and exit 


of the test section, local friction factors could not be determined in this 


i | 
investigation. Overall average friction factors with heat addition were 
compared to the correlation proposed by Taylor [197] 
= -0.32 -0.5 
Fay = (0.0014 + 0.125Re,, (Ty, av/Th av? (10) 
| 
| 
for the data of a wide variety of experiments with gas flow with 
0.62 < Pr < 0.81. Most previous measurements agreed within 10%. For 
,° 


31 


evaluation of this expression, integrated averages of both the local wall 
and local bulk gas temperature were used along with the average pressure to 
determine average density and viscosity. The overall friction factor was 


determined from the frictional pressure drop, 


T T 
A os = _ GPR es ee (11) 
Pey ~ Py ~ Po 5. Py 


and the modified wall Reynolds number was defined as 


Re, mS (GD/H) (Ty ay/Twyav? (12) 


In the present study correlation (10) predicted most of the data within 52. 
These measurements are presented in Figure 7. The few data deviating more 
than -5% from the correlation were mixture runs with maximum wall 
temperatures in excess of 668°C (1234°F) and T,/T, > 1.8 and with modified 
wall Reynolds numbers in the low range between 8,000 and 20,000. The data 
deviating more than +5Z% were all air runs with modified wall Reynolds 
numbers between 8,500 and 17,500. There appeared to be no effect of 
Prandtl number. 

Experimental local Nusselt numbers were compared to Nusselt numbers 
predicted by modified forms of the prediction equations recommended by 


Pickett, Taylor and McEligot [1979] 


0.8 0.55 0.4 
= ; 3 
Nu, = 0.021 Re," ” Pry (T/T) + 0.85 D/x] (13) 
and by Taylor [1%8] 
= 0.8 0.4 ~a (14) 
Nu, = 0.023 Re. Pry (T,,/T,) 
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where a= (0.57 - — 


x/D | 


For the present comparisons a value of 0.65 was finally taken as the 
exponent of the Prandtl number in both of these equations. 

Figure 8 shows the ratios of the measured Nusselt number to 
the Nusselt number predicted from Equation (13) and (14) as functions 
of x/D for two heating rates of air, He - Xe (M = 14.5, Pr = 0.30) and 
Hy - Xe(M = 29, Pr = 0.18). The data shown here are typical of all the 
data measured in this investigation. Both of the correlation equations 


appear to predict downstream Nusselt numbers with acceptable accuracy for 


| 
| 
| 
| 
| 
all three Prandtl numbers. Both equations predict entrance effects better 
for the low heating rate than for the higher rate and are in close 

agreement with each other, even though the methods of handling entrance 

effects are different. A better understanding of entrance effects is 

needed. 

Increasing the exponent of the Prandtl number to 0.65 yielded 
considerable improvement and promise in the predictions using Equations (13) | 
and (14) for all Prandtl numbers. The value of Pr°65/pr9-4 is 0.65 for 
Pr = 0.181, 0.74 for Pr = 0.30 and 0.92 for air; it is obvious that | 
determining the correct exponent for Prandtl number is far more difficult | 
for values of Prandtl number near 1 than in the lower Prandtl number cases, 
where the sensitivity of the Nusselt number to the Prandtl number exponent 


is much greater. 
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Examination of air data with Reynolds numbers of ~ 34,000, 65,000 
and 85,000 showed less than 3% variation of Nu/Nuy with Reynolds number for 
a particular heating rate for all x/D greater than 2.2. The same held true 
for the gas mixtures when two Reynolds numbers were available for 


comparison at the same mixture concentration. 


CONCLUSION 


From experiments with heated flow of helium-xenon and hydrogen- 
Xenon mixtures with Prandtl numbers in the range 0.18 to 0.30, the 
following major conclusions may be drawn. The Colburn analogy and the 
Dittus-Boelter correlation based on measurements for Pr > C.7, greatly 
overpredict Nusselt numbers for fully established conditions with constant 
properties in this range. The correlation equation of Sleicher and Rouse 
and also that of Churchill underpredict the Nusselt number in this range. 
Of the correlations examined, that of Petukhov [1970] best represents the 
data the constant property Nusselt number at 0.18 < Pr< 0.72 in fully 
established conditions. 

The friction coefficients for flow with property variation 
were predicted within ts¢ by the correlation proposed by Taylor [1967], 
Equation (10), for modified wall Reynolds numbers greater than 20,000. 

Nusselt numbers for flow with property variation were predicted 
using the correlation recommended by Pickett, Taylor and McEligot [1979], 
Equation (13), and Taylor [1968], Equation (14), but with the exponent of 


Prandtl number increased 0.65 in both cases. The predictions of both 
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4 


correlations are in good agreement with each other. For low heating rates 
the measured and predicted values agreed within “10% for axial distances 
from 4.3 to 52.3 diameters. For the higher heating rates both equations 
failed to predict the effects of the entrance region adequately at x/D less 
than about 20. Most of the data at x/D greater than 20 was predicted 


within +102. 
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APPENDIX A 
Uncertainty Analysis 


To determine the validity of the deduced results it is necessary to 
determine their estimated experimental uncertainties. An analysis has been 
conducted in accordance with recognized procedures to estimate the percent 
uncertainty in the important heat transfer and friction characteristics 
along the test section. Doebelin [1%6] considers the problem of computing 


a quantity N, where N is a known function of the n independent variables, 


U1, u2; U3; eee y Un. That is, 


N = fn (uy, Ups bas ee es u, (Al) 


iC The u's are the measured quantities (instrument or component outputs) and 
are uncertain by +Au), +Aug, +\u3, ... , +\Un, respectively. These 
uncertainties will cause an uncertainty AN in the comput2d result N. We 

Cd & are concerned with experimental uncertainties here rather than systematic 
errors, like the thermocouple conduction error which can be corrected. The 
u's may be considered as absolute limits on the uncertainties, as 
statistical bounds such as 39 limits, or as uncertainties on which we are 
willing to give certain odds as including the actual error. However, the 
method of computing AN and the interpretation of its meaning are different 


for the first case as compared with the second and third. After 


calculating the limits on the uncertainty in N, then N + AN and the 


x 
» 
» 
. 
ie 


wre 


percentage uncertainly is known. In all cases, systematic errors (bias) 
were removed by calibration where they were known to exist. 
If the individual errors are thought of as +30 limits, then, as 


shown by Kline and McClintock [1953], the general equation used is 


AN =E = Hou 2 + (Au ot )2 +... + (Au euleye (A2) 
*y 


arss 2du nou 
2 n 


where E..g3 represents a +30 limit on N and 99.7 percent of the values of N 
can be expected to fall within these limits (arss = absolute root-sum 
square). 

An analysis was performed to determine the uncertainty of local 
bulk Nusselt number, Nu, as calculated from the measured experimental data. 
Table Al lists the uncertainties of the instruments (Au's) used in this 
investigation. The uncertainties of the directly measured quantities were 
determined from manufacturers’ specifications and experience, The 
uncertainties of the gas properties were not included. 

The estimated uncertainty in the experimental Nusselt number which 
results from the uncertainty in a measured variable {is shown in Table A2. 
The measured variables included mass flow rate, current, inlet bulk 
temperature of the gas, wall temperature, inlet static pressure, and the 
resistance per unit length of the test section. Uncertainties, Au, in 
these measurements were assigned from Table Al and their values are shown 
in Table A2. An uncertaicty in mass flow rate of 1.5 percent was used in 
Table A2. The values of 3} Nu/ju in Table A2 were attained by changing the 


variable by a: amount equal to its estimated uncertainty and rerunning the 
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computer program for data reduction at the experimental conditions. It was 


found that a change of 0.2 psia to the inlet static pressure did not change 
the experimental Nusselt number noticeably. Also, an uncertainty of 
0.00008 Q/in in the measurement of resistance per unit length had an 
insignificant effect on the percentage uncertainty in the experimental 
Nusselt number. 

Three types of comparisons can be made from Table A2. A comparison 
between the two runs shows that the percentage uncertainty for a relatively 
low heating rate run to be higher than a high heating rate run at about the 
same Reynolds number. The result is due to the relatively high percentage 
uncertainties in the inlet bulk and wall temperatures being much more 
dominant in the lower heating rate case. Secondly, the change in the 
percentage uncertainty along the axial length of the test section can be 
examined. For the low heating rates the major uncertainties in the Nusselt 
number are provided by the inlet bulk gas temperature and the wall 
temperatures; they decrease with x/D. For higher heating rates, the 
uncertainty increases with x/D because the uncertainties in tube wall 
temperature and, therefore, temperature difference increase significantly 
with temperature level, while the contributions of uncertainties in mass 
flow rate, electrical power and inlet gas temperature remain small. 
Thirdly, the individual variables can be compared to each other for the two 


different heating rates as well as for the range of x/D. 


APPERDIX B 


EXPERIMENTAL DATA WITH HEAT ADDITION 


Table Bl. Summary of Experimental Data 


Run Gas " Pri Re; (Ty/Tp) max q*max 
686H He-Xe 83.8 0.25 85147 1.18 0.0007 
687H 85802 1.42 0016 
6888 86740 1.65 . 0026 
68 9H 87373 1.89 .0035 
6 90H 86063 2.17 0049 
691H 32646 1.16 .0009 
6 92H 32647 1.39 0022 
6 3H 33513 1.62 0034 
6 94H 33328 1.88 0051 
695H 33947 2.22 0069 
6 96H He-Xe 40.0 0.21 36 183 1.99 .0053 
6 97H 34994 1.75 0041 
6 98H 34280 1.59 0032 
699H 34686 1.36 .0018 
7008 34865 1.15 0007 
701H 58657 1.83 .0039 
702H 58525 1.59 .0027 
703H 61534 1.38 .0017 
7048 59851 1.17 0008 
705H 61772 1.77 .0037 
707H He-Xe 28.3 0.23 53390 1.67 0030 
708H 55449 1.38 0017 
7098 48411 2.06 0051 
710H 49001 1.16 0007 
7118 He-Xe 14.5 0.30 40869 1.15 0006 
7128 40554 1.32 0013 
713H 35374 1.61 0027 
7148 38 962 1.84 0034 
715H 34042 2.04 .0047 
7198 Hy-Xe 29.0 0.18 73916 1.48 0.0021 
7208 71174 1.78 0034 
7218 72531 1.30 0013 
722H 73803 1.13 0005 


Table B1. Summary of Experimental Data 


t (continued ) 

C 

Re Run Gas M Pr; Re; (Ty /Th) max q*+max 
72348 Air 28.97 0.72 84169 1.20 -0005 
7248 84455 1.19 -0005 
7258 84778 1.46 0012 
7268 64990 1.40 -0012 
727H 66076 1.19 0005 
7288 34149 1.74 -0023 
7298 34946 1.45 -0014 
7308 34765 1.20 -0006 
7318 33882 2.03 -0033 
7328 34777 2.35 0044 
7334 84294 1.74 0019 
7348 84742 2.00 -0026 
7358 85824 2.30 -0033 
736H 64311 1.90 -0024 
73748 66215 1.73 0019 
7388 64199 2.38 -0038 


The headings and their definitions used in the listing of the heated flow 


data are below. 


Heading Definition 
TIN Inlet gas temperature, °F 
TOUT Calculated outlet gas temperature, °F ; 
I Alternating current, amperes 
E Voltage drop between voltage taps, volts 
PR, IN Inlet Prandtl number, ¢, Wik ji 
GR/RESQ Ratio of Grashof number g D4q" v7 Kr), to the 

square of the inlet Reynolds number, GD/ Uj 
MACH (2) Mach number at thermocouple 2 
MACH (16) Mach number at thermocouple 16 


Me tat 2 = oa on —— 


Heading 
TSURR 
Q+(8) 


TC 


HICOEF 


BULK REYNOLDS 


BULK NUSSELT 
PT 

TB 

PRESS DEFECT 
TW,AV(F) 
TB,AV(F) 
TW,AV/TB,AV 
DELTA P(PSI) 


RE,B,AV 


RE,W,MOD,AV 


F,B,AV 


Definition 
Temperature of surroundings inside draft shield, °F 
Nondimensional turbulent heat flux parameter. Corresponds 
to q* in text at thermocouple 8, qt = q"/(G Cpt T;) 
Thermocouple number 
Axial position, corresponds to x/D in text 
Ratio of heat loss to heat flux to gas 
Inside tube wall temperature, °F 
Wall to bulk temperature ratio 
Heat flux to gas, Btu/hr-ft2 
Heat transfer coefficient, Btu/hr-ft2-°F 
Reynolds number evaluated at bulk temperature, GD/yy, 
Nusselt number evaluated at bulk temperature, hD/kp 
Pressure tap: l-near inlet, 2-near outlet 
Bulk static temperature °F 
Pressure defect, 0; g. (p; - p)/G? 
Average inside tube wall temperature °F 
Average bulk gas temperature °F 
Average wall to average bulk temperature ratio 
Pressure drop from start of heating to pressure tap #2 


Average value of bulk Reynolds number from start of 
heating to pressure tap #2, GD/u, ae 
tT 


Average value of modified wall Reynolds number, 
GD/Hy ay (Ty av/Ty, ay) 


Friction factor through heated section of tube, 
BcPb ay D Pe, /2LG 


A A ES 
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